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Abstract. This is a review on cosmological perturbation theory. After an intro- 
duction, it presents the problem of gauge transformation. Gauge invariant variables 
are introduced and the Einstein and conservation equations are written in terms of 
these variables. Some examples, especially perfect fluids and scalar fields are pre- 
sented in detail. The generation of perturbations during inflation is studied. Lightlike 
geodesies and their relevance for CMB anisotropics are briefly discussed. Perturba- 
tion theory in braneworlds is also introduced. 
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1 Introduction 



The idea that the large scale structure of our Universe might have grown our 
of small initial fluctuations via gravitational instability goes back to Newton 
(letter to Bcntlcy, 1692[1]). 

The first relativistic treatment of linear perturbations in a Friedmann- 
Lemaitre universe was given by Lifshitz (1946) [2]. There He found that the 
gravitational potential cannot grow within linear perturbation theory and he 
concluded that galaxies have not formed by gravitational instability. 

Today we know that it is sufficient that matter density fluctuations can 
grow. Nevertheless, considerable initial fluctuations with amplitudes of the 
order of 10~^ are needed in order to reproduce the cosmic structures observed 
today. These are much larger than typical statistical fluctuations on scales of 
galaxies and we have to propose a mechanism to generate them. Furthermore, 
the measurements of anisotropies in the cosmic microwave background show 
that the amplitude of fluctuations is constant over a wide range of scales, the 
spectrum is scale independent. 

As we shall see, standard inflation generically produces such a spectrum 
of nearly scale invariant fluctuations. 

In this course I present gauge invariant cosmological perturbation theory. 
I shall start by defining gauge invariant perturbation variables. Then I present 
the basic perturbation equations. As examples for the matter equations we 
shall consider perfect fluids and scalar fields. The we briefly discuss lightlike 
geodesies and CMB anisotropies (this section will be very brief since it is 
complemented by the course on CMB anisotropies by A. Challinor). Finally, 

1 shall make some brief comments on perturbation theory for braneworlds, a 
topic which is still wide open in my opinion. 

2 The background 

I shall not come back to the homogeneous universe which as been discussed 
in depth in the course by K. Tamvakis. I just specify our notation which is as 
follows: 

A Friedmann-Lemaitre universe is a homogeneous and isotropic solution 
of Einstein's equations. The hyper-surfaces of constant time are homogeneous 
and isotropic, i.e. , spaces of constant curvature with metric {r])jij dx^ dx^ , 
where 7^ is the metric of a space with constant curvature k. This metric can 
be expressed in the form 




dr^ + X^ir) {d-d^ + sin^Mip^) 
( r , K = 



(1) 
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where we have rescaled a{r]) such that k = ±1 or 0. (With this normahzation 

the scale factor a has the dimension of a length and rj and r are dimensionless 
for K ^ 0.) The four-dimensional metric is then of the form 

gn^cix'^dx" = -a^{ri)drf' + a^{ri)^ijdx'dx' . (3) 

Here rj is called the conformal time. The physical or cosmological time is given 
by dt = adt]. 

Einstein's equations reduce to ordinary differential equations for the func- 
tion a{ni) (with ' = d/drj): 



+ K = n'^ + K= ^^a^p + ^Aa^ (4) 



a \ 47rG 2/ o\ ^ a I 0'\ I cl 



2 



where p = — Tg , p = T? (no sum!) and all other components of the energy 
momentum tensor have to vanish by the requirement of isotropy and homo- 
geneity. A is the cosmological constant. We have introduced H = d/a. The 
Hubble parameter is defined by 

a a-' 

Energy momentum "conservation" (which is also a consequence of (4) and 
(5) due to the contracted Bianchi identity) reads 

p = -3(^^yp + p) = -3{l + w)n, (6) 

where w = p/p. Later we will also use = p/ p. From the definition of w and 
p together with Eq. (6) one finds 

w = 3{1 + w){w ~ cl)n . (7) 

From the Friedmann equations one easily concludes that for k = yl = and 
w =const. the scale factor behaves like a power law, 

2 2 

a oc 77i+3» oc i^cn-"-) . (8) 

Important examples are 

q = 2 for dust, iv = 

accrf with {(7=1 for radiation, w = 1/3 (9) 

q = —1 for inflation (or a cosm. const.), w = —1 



We also define 
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3^2 
Q = — 



(10) 



where the index o indicates the value of a given variable today. Friedmann's 
equation (4) then requires 



One often also uses 12 = i?^ + Ha = ! — 



(11) 



3 Gauge invariant perturbation variables 

The observed Universe is not perfectly homogeneous and isotropic. Matter is 
arranged in galaxies and clusters of galaxies and there are large voids in the 
distribution of galaxies. Let us assume, however, that these inhomogeneities 
lead only to small variations of the geometry which we shall treat in first order 
perturbation theory. For this we define the perturbed geometry by 

9t,u = gij.v + ea^h^^ (12) 

Qfj^i, being the unperturbed Priedmann metric. We conventionally set (absorb- 
ing the "smallness" parameter e into /i^,^) 

= Sniy + a^h^v, goo = -a^, 9ij = a^Hj \hfi„\^l 
_ _„ . (13) 

Tjf = T';:+ei^, To = -p, Tj=psj m/p^i. 



3.1 Gauge transformation, gauge invariance 

The first fundamental problem we want to discuss is the choice of gauge in 
cosmological perturbation theory: 

For linear perturbation theory to apply, the spacetime manifold A4 with 
metric g and the energy momentum tensor T of the real, observable universe 
must be in some sense close to a Friedmann universe, i.e. , the manifold M 
with a Robertson-Walker metric g and a homogeneous and isotropic energy 
momentum tensor T. It is an interesting, non trivial unsolved problem how 
to construct 'the best' g and T from the physical fields g and T in practice. 
There are two main difficulties: First, spatial averaging procedures depend 
on the choice of a hyper-surface of constant time and they do not commute 
with derivatives, so that averaged fields g and T will in general not satisfy 
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Einstein's equations. Secondly, averaging is in practice impossible over super- 
horizon scales. 

Even though we cannot give a constructive prescription, we now as- 
sume that there exists an averaging procedure which leads to a Priedmann 
universe with spatially averaged tensor fields Q, such that the deviations 
{Tfj,u-T / va.&-K{af3}{\T anW and (3^1. -5^^,)/ max{„^}{5„^} arc small, and 5 
and T satisfy Priedmann's equations. Let us call such an averaging procedure 
'admissible'. There may be many different admissible averaging procedures 
(e.g. over different hyper-surfaces) leading to slightly different Priedmann 
backgrounds. But since \g — g\ is small of order e, the difference of the two 
Priedmann backgrounds must also be small of order e and we can regard it as 
part of the perturbation. 

We consider now a fixed admissible Priedmann background {g, T) as cho- 
sen. Since the theory is invariant under diffeomorphisms (coordinate transfor- 
mations), the perturbations are not unique. Por an arbitrary diffeomorphism 
(j) and its push forward </>*, the two metrics g and (lf{g) describe the same 
geometry. Since we have chosen the background metric g we only allow dif- 
feomorphisms which leave g invariant i.e. which deviate only in first order 
form the identity. Such an 'infinitesimal' diffeomorphism can be represented 
as the infinitesimal flow of a vector field X , cf) ^ (fif. Remember the defini- 
tion of the flow: Por the integral curve 7a; (s) of X with starting point x, i.e., 
7x(s = 0) = a; we have (/>^(x) — -fxis)- In terms of the vector field X, to first 
order in e, its puUback is then of the form 

(f)* = id + eLx 

{Lx denotes the Lie derivative in direction X). The transformation g (j)*{9) 
is equivalent to g+ea^h g+e{a'^h+Lxg)+0{€^), i.e. under an 'infinitesimal 
coordinate transformation' the metric perturbation h transforms as 

h^h + a-'^Lxg. (14) 

In the context of cosmological perturbation theory, infinitesimal coordinate 

transformations are called 'gauge transformations'. The perturbation of a ar- 
bitrary tensor field Q = Q + eQ'^^^ obeys the gauge transformation law 

qW^qW+LxQ. (15) 

Since every vector field X generates a gauge transformation (f) = (jy^ , we 
can conclude that only perturbations of tensor fields with LxQ = for all 
vector fields X, i.e., with vanishing (or constant) 'background contribution' 
arc gauge invariant. This simple result is sometimes referred to as the 'Stewart- 
Walker Lemma' [3]. 

The gauge dependence of perturbations has caused many controversies in 
the literature, since it is often difficult to extract the physical meaning of gauge 
dependent perturbations, especially on super-horizon scales. This has led to 



6 Ruth Durrer 



the development of gauge invariant perturbation theory which we are going to 

use throughout this review. The advantage of the gauge -invariant formalism is 
that the variables used have simple geometric and physical meanings and are 
not plagued by gauge modes. Although the derivation requires somewhat more 
work, the final system of perturbation equations is usually simple and well 
suited for numerical treatment. We shall also see, that on sub-horizon scales, 
the gauge invariant matter perturbation variables approach the usual, gauge 
dependent ones. Since one of the gauge invariant geometrical perturbation 
variables corresponds to the Newtonian potential, the Newtonian limit can be 
performed easily. 

First wc note that since all rclativistic equations arc covariant (i.e. can be 
written in the form Q = for some tensor field Q), it is always possible to 
express the corresponding perturbation equations in terms of gauge invariant 
variables [4, 5, 6]. 

3.2 Harmonic decomposition of perturbation variables 

Since the {77 = const} hyper-surfaces are homogeneous and isotropic, it is sen- 
sible to perform a harmonic analysis: A (spatial) tensor field Q on these hyper- 
surfaces can be decomposed into components which transform irreducibly un- 
der translations and rotations. All such components evolve independently. 
For a scalar quantity / in the case k = this is nothing else than its Fourier 
decomposition: 

/(x,^?) = ydV>)e*^ (16) 

(The exponentials lk;(x) = e^^^ are the unitary irreducible representations of 
the Euclidean translation group.) For k = 1 such a decomposition also exists, 
but the values k are discrete, fc^ = l{i-\-2) and for k = —1, they are bounded 
from below, k"^ > 1. Of course, the functions Yk are different for k^O. 

They form the complete orthogonal set of eigenfunctions of the Laplacian, 

AY^^^ = -k'^Y^^\ (17) 

In addition, a tensorial variable (at fixed position x) can be decomposed 
into irreducible components under the rotation group SO{3). 

For a vector field, this is its decomposition into a gradient and a rotation, 

Vi = Viif + Bi, (18) 

where 

B\i = 0, (19) 

where we used X|j to denote the three-dimensional covariant derivative of X. 
Here is the spin and B is the spin 1 component of the vector field V. 
For a symmetric tensor field we have 
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= JfL7i,- + (v. V, - izi7^i) Ht + \ (ifjy^ + H^^^) + h\J^ , (20) 



^v),.^^Tr^^T,^0_ (21) 



where 

Here iJi and Ht are spin components, H^^^ is a spin 1 component and 

(T) 

H^j is a spin 2 component. 

Wc shall not need higher tensors (or spinors). As a basis for vector and 
tensor modes we use the vector and tensor type eigenfunctions of the La- 
placian, 

= -k^yp and (22) 
AyIP = -k'YjP , (23) 

where Y^^^ is a transverse vector, Yj^ =0 and Yjp is a symmetric trans- 
verse traceless tensor, Yj-'^^-' — YjP^^ = 0. 

According to Eqs. (18) and (20) we can construct scalar type vectors and 
tensors and vector type tensors. To this goal we define 

Y.(^) = -fc-iyj^.^) (24) 
Y!f^ ^ fc-l^(f + Ij,,Y^s^ (25) 

In the following wc shall extensively use this decomposition and write down 
the perturbation equations for a given mode k. 

The decomposition of a vector field is then of the form 

B, = ByI^'> + b'^^^Y^^K (27) 

The decomposition of a tensor field is given by (compare Eq. (20)) 

H,j = H^Y^^hij + HtY^P + H(^)yP + if (^)Ff ). (28) 

Here B, S(^), Hl, Ht, H^^^ and if(^) are functions of r? and k 

3.3 Metric perturbations 

Perturbations of the metric are of the form 

g^iu = 9p,u + a^h^^. (29) 

We parameterize them as 
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h^^^dxi^dx" = -2Adrf - 2Bidr]dx' + 2Hijdx'dx^, (30) 

and we decompose the perturbation variables Bi and Hij according to (27) 
and (28). 

Let us consider the behavior of ft.^^ under gauge transformations. We set 
the vector field defining the gauge transformation to 

X=^Tdn + L^di. (31) 

Using simple identities from differential geometry like L-x^{df) = d{L-x^f) and 
{Lx.l)i3 = ^i\3 + we obtain 



ix5 = a 



2 i^T + t\ drf- +2 {Li - Ti) dr}dx' 



+ {^'^T-iij + Li\j + Lj\^ dx'dx^ 



(32) 



Comparing this with (30) and using (14) we obtain the following behavior 
of our perturbation variables under gauge transformations (we decompose the 
vector Li = Ly/^)+L(^)y/^)): 

A+-T + f (33) 
a 

B ^ B-L-kT (34) 
Biy) ^ B(^) - L(^) (35) 

Hl^Hl + -T+^L (36) 
a 6 

Ht->Ht- kL (37) 

H(y) ^ F(^) - fci(^) (38) 

Two scalar and one vector variable can be set to zero by a cleverly chosen 
gauge transformations. 

One often chooses kL = Ht and kT = B + L, so that the variables Ht 
and B vanish. In this gauge (longitudinal gauge), scalar perturbations of the 
metric are of the form {Hx = B = 0): 

h^s) = _2a^dr,2 _ 2$^ijdx'dxK (40) 

'1' and # are the so called Bardeen potentials. In a generic gauge they are 
given by 

^ = A--k-'^a-k-^& (41) 

a 

$ = -Hl - \ht + -k-^a (42) 
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with a = k~^HT — B. A short calculation using Eqs. (33) to (37) shows that 
!Z' and <P arc gauge invariant. 

In a Friedmann universe the Weyl tensor vanishes. It therefore is a gauge 
invariant perturbation. For scalar perturbations one finds 



(43) 



All other components vanish. 

For vector perturbations it is convenient to set kL^^^ = H^^^ so that H^^^ 
vanishes and we have 

h^Ux^'dx" = 2a^^'^Yl^^dr]dx\ (44) 

We shall call this gauge the "vector gauge". In general cr^^^ = k~^H^^^ —B^^^ 
is gauge invariant^. The Weyl tensor from vector perturbation is given by 

C^r. = l-'^'HYj^i - - ^7..>112'= + p^mYir'] (46) 

Clearly there are no tensorial (spin 2) gauge transformation and hence H^J^ is 
gauge invariant. The expression for the Weyl tensor from tensor perturbation 
is identical to the one for vector perturbation upon replacement of cr(J^'' by 

3.4 Perturbations of the energy momentum tensor 

Let T,^ = t'^ + 9^^ be the full energy momentum tensor. We define its en- 
ergy density p and its energy flow 4-vector u as the time-like eigenvalue and 
eigenvector of T^f: 

Tl^u" = -puf", V? = -1. (47) 
We then parameterize their perturbations by 

p = p(l + ,5), M = M°5t+M*5,. (48) 

u° is fixed by the normalization condition, 

u^ = \\-A). (49) 
a 

We further set 

= -v' = ?;y(^)' + (50) 
a 

^ Yij^ ^a''^-' is the shear of the hyper-surfaces of constant time. 
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We define = u^^Ui, + 6^, the projection tensor onto the part of tangent 
space normal to u and set the stress tensor 

T^"" = P^P^T"^. (51) 

In the unperturbed case we have Tq = 0, rj = p6j. Including perturbations, 
to first order we still obtain 

=r°=T^ = 0. (52) 
But rj contains in general perturbations. We set 

rj = p [(1 + ttl) 5} + 77j] , with ni = 0. (53) 
We decompose IJj as 

TTj = n^s)yis) ^ ^ jj(v)Yiy) ^ + n^T)YiT)i_ (54) 

We shall not derive the gauge transformation properties of these pertur- 
bation variables in detail, but just state some results which can be obtained 
as an exercise (see also [5]): 

• Of the variables defined above only the n'^"^'^''^^ are gauge invariant; they 
describe the anisotropic stress tensor, IT^ = — ^st^S^. They are gauge 
invariant due to the Stewart-Walker lemma, since H = 0. For perfect 

fluids nif = 0. 

• A second gauge invariant variable is 

r = TTi - ^S, (55) 

w 

where = p/p is the adiabatic sound speed and w =p/p is the enthalpy. 
One can show that F is proportional to the divergence of the entropy flux 
of the perturbations. Adiabatic perturbations are characterized by r" = 0. 

• Gauge invariant density and velocity perturbations can be found by com- 
bining S, V and v^^^ with metric perturbations. 



We shall use 



V = v- \ht = v''^°''s) (56) 
k 

D, = 6 + 3(1 + w)^{k-^HT - k-^B) = (5(i°°s) (57) 

£) = ^(lons)+3(l+^) ("^^ (58) 

\a J k 

Dg = 5 + 3(1 + w) (hl + ^Ht\ = _ 3(1 + (59) 
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V^y) = - lif(^) = ^(vec) (60) 

Q = - B(^) = - B(^) (61) 

Here ti(^°"s), ijCiong) and w^™'^'' are the velocity (and density) perturbations in 
the longitudinal and vector gauge respectively, and cr^^' is the metric per- 
turbation in vector gauge (see Eq. (44)). These variables can be interpreted 
nicely in terms of gradients of the energy density and the shear and vorticity 
of the velocity field [7] . 

Here I just want to show that on scales much smaller than the Hubble scale, 
krj ^ 1, the metric perturbations are much smaller than 6 and v and we can 
thus "forget them" (which will be important when comparing experimental 
results with calculations in this formalism): 

The perturbations of the Einstein tensor are given by second derivatives 
of the metric perturbations. Einstein's equations yield the following order of 
magnitude estimate: 

= o(^-^h+^h + k^hj (63) 




= 0{h + kr]h+ (krjfh) . (64) 

For kr}^ 1 this gives 0{5,v) = O (^) ^ 0{h). On sub-horizon scales the 
difference between 5, and D is negligible as well as the difference 

between v and V or v^^', y(^) and 



4 Einstein's equations 

We do not derive the first order perturbations of Einstein's equations. This 
can be done by different methods, for example with Mathematica. We just 
write down those equations which we shall need later. 



4.1 Constraint equations 



4nGa'^pD = -{k^-3K)^ (00) | 
4nGa\p + p)V = k ((f) ^ + S) (Oi) | ^''^''^^'^ ^^^^ 

SwGa^ip + p)n=^ {2k - fc^) CT^^' (Oi) (vector) (66) 
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4.2 Dynamical equations 



SirGa^pn'^^'i (scalar) (67) 
87rGa2p77(^) (vector) (68) 

SirGa^pnlp (tensor) (69) 

There is a second dynamical scalar equation, which is somewhat combersome 

and not really needed, since we may use one of the conservation eqs. given 
below instead. Note that for perfect fluids, where TIj = 0, we have = 
£7^^' a 1 /a^ and H^"^^ obeys a damped wave equation. The damping term can 
be neglected on small scales (over short time periods) when r]~^ ~ 2k + k^, 
and Hij represents propagating gravitational waves. For vanishing curvature, 

these are just the sub-horizon scales, /cry ~ 1. For k <Q, waves oscillate with 
a somewhat smaller frequency, u) = \J2k + fc^, while for /t > the frequency 
is somewhat larger than fc. 



A;2 {^-^) = 



4.3 Energy momentum conservation 

The conservation equations, T.^" = lead to the following perturbation equa- 
tions. 



► (scalar) (70) 



Dg + 3 (c2 - w) (i) Dg + {1 + w)kV + 3w (|) r = 
V + (t) (1 - 3c^) V = k{>F + 3cl<P) + ^D, 
+^ [^-§(1-1?)^] 

+ (1 - 3c^) (^) = ^ - I ) nr (vector) (71) 

These can of course also be obtained from the Einstein equations since they 
are equivalent to the contracted Bianchi identities. For scalar perturbations 
we have 4 independent equations and 6 variables. For vector perturbations 
we have 2 equations and 3 variables, while for tensor perturbations we have 
1 equations and 2 variables. To close the system we must add some matter 
equations. The simplest prescription is to set F = 77jj = 0. This matter 
equation, which describes adiabatic perturbations of a perfect fluid gives us 
exactly two additional equations for scalar perturbations and one each for 
vector and tensor perturbations. 

Another simple example is a universe with matter content given by a 
scalar field. We shall discuss this case in the next section. More complicated 
examples are those of several interacting particle species of which some have 
to be described by a Boltzmann equation. This is the actual universe at late 
times, say z ~ 10''. 
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4.4 A special case 

Here we want to rewrite the scalar perturbation equations for a simple but 
important special case. We consider adiabatic perturbations of a perfect fluid. 
In this case 77 = since there are no anisotropic stresses and 7^ = 0. Eq. (67) 
then implies <!> = ^. Using the first equation of (65) and Eqs. (59,58) to replace 
Dg in the second of Eqs. (70) by ^ and V, flnally replacing V by (65) one can 
derive a second order equation for 'P, which is, in this case the only dynamical 
degree of freedom 

If- + 3H(1 + cl)^ + [(1 + 3c^)(H2 - k) - (1 + 2,w)n^ + clk'^]^ = . (72) 

Another interesting case (especially when discussing inflation) is the scalar 
field case. There, as we shall see in Section 6, iT = 0, but in general F ^ 
since Sp/dp ^ p/p. Nevertheless, since this case again has only one dynamical 
degree of freedom, wc can express the perturbation equations in terms of one 
single second order equation for ]P. In Section 6 we shall find the following 
equation for a perturbed scalar field cosmology 

W + 3n{l + cl)>P + [(1 + 3c2)(h2 - k) - (1 + 3w)n^ + fc2]if = . (73) 

The only difference between the perfect fiuid and scalar field perturbation 
equation is that the latter is missing the factor in front of the oscillatory 

term. Note also that for k = and w = = constant the time dependent 
mass term m'^{r]) = — (1 + 3cl){Tl'^ — k) + {1 + 3w)Ti? vanishes. It is useful to 
define also the variable [11] 

r n-i/2 
u = a AnGiH'^ - H + k) ^, (74) 

which satisfies the equation 

it + {Tk^ ~ e/e)u = 0, (75) 

where T = or T = 1 for a perfect fluid or a scalar fleld background respec- 
tively, and 

0= • (76) 

Another interesting variable is 

Using Eqs. (72) and (73) respectively one obtains 



(78) 
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hence on super horizon scales, k/7i <C 1, this variable is conserved. 

The evolution of C is closely related to the canonical variable v defined by 

. = -"^^C. (79) 

which satisfies the equation 

V + {Tk^ - z/z)v = , (80) 

for 

TH • (^'^ 

More details on the significance of the canonical variable v will be found in 
sections 6 and 7. 



5 Simple examples 

We first discuss two simple applications which are important to understand 
the CMB anisotropy spectrum. 

5.1 The pure dust fluid for k = 0,A = 

We assume the dust to have w = = p = and 11 = F = 0. Equation (72) 
then reduces to 

^+-<ir = Q (82) 

V 

with the general solution 

<Z' = >fo + >fi4 (83) 

with arbitrary constants i'o and •Z'l . Since the perturbations are supposed to be 
small initially, they cannot diverge for r] ^ 0, and we have therefore to choose 
the growing mode, = 0. Another way to argue is as follows: If the mode 
tf'i has to be small already at some early initial time ryin, it will be even much 
smaller later and may hence be neglected at late times. But also the if'o mode is 
only constant. This fact led Lifshitz who was the first to analyze cosmological 
perturbations to the conclusions that linear perturbations do not grow in a 
Friedman universe and cosmic structure cannot have evolved by gravitational 
instability [2]. However, the important point to note here is that, even if the 
gravitational potential remains constant, matter density fluctuations do grow 
on sub-horizon scales and therefore inhomogeneities, structure can evolve on 
scales which are smaller than the Hubble scale. To see that we consider the 
conservation equations (70), (67) and the Poisson equation (65). For the pure 
dust case, w = cl = n = r = 0, they reduce to 
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bg = —kV (energy conservation) 
V + (^-^V = k\[' (gravitational acceleration) 

-k'^^ = AnGa'^p (^Dg + 3 + fc-^y 

where we have used the relation 



D = Dg + 3{l + w)[-^+{-] k-'V 



a. 

The Friedmann equation for dust gives AirGpa"^ = 3/2(^/0)^ 
kr] = X and ' = d/dx, the system (84-86) becomes 



D'g = -V 

X 



L»o + 3( <f + -1/ 



-If. 



(84) 
(85) 

(Poisson), (86) 



(87) 

6/77^. Setting 

(88) 
(89) 

(90) 



We use (90) to eliminate ^ and (88) to eliminate Dg, leading to 

(18 + y. + + v'-(^-^^+Ay = Q. (91) 



The general solution of Eq. (91) is 



1/ = Fox + 4 

x^ 



(92) 



The Vi mode is the decaying mode (corresponding to ^'i) which we neglect. 
The perturbation variables are then given by 



V = Vox 



Wc distinguish two regimes: 
i) super-horizon, x <C 1 where we have 







V = 


l^x 


Dg = 




^ = 


% . 



(93) 
(94) 
(95) 



(96) 

(97) 
(98) 



Note that even though V is growing, it always remains much small than ^ 
or Dg on super-horizon scales. Hence the largest fluctuations are of order <]/ 
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which is constant. 

ii) Sub-horizon, a; 1 where the solution is dominated by the terms 

V = (99) 

Dg = -l^ox'' (100) 
If- = !Z>o = constant. (101) 

Note that for dust 

D = Dg + M+-V = -l%x'' . 

X 6 

In the variable D the constant term has disappeared and we have D <i; •Z' on 
super-horizon scales, a; ^ 1 . 

On sub-horizon scales, the density fluctuations grow Hke the scale factor 
<x x^ <x a. Nevertheless, Lifshitz' conclusion [2] that pure gravitational in- 
stabiHty cannot be the cause for structure formation has some truth: If we 
start from tiny thermal fluctuations of the order of 10~^^, they can only grow 
to about 10~^° due to this mild, power law instability during the matter 
dominated regime. Or, to put it differently, if we want to form structure by 
gravitational instability, we need initial fluctuations of the order of at least 
10~^, much larger than thermal fluctuations. One possibility to create such 
fluctuations is quantum particle production in the classical gravitational field 
during inflation. The rapid expansion of the universe during inflation quickly 
expands microscopic scales at which quantum fluctuations are important to 
cosmological scales where these fluctuations are then "frozen in" as classical 
perturbations in the energy density and the geometry. We will discuss the 
induced spectrum on fluctuations in Section 7. 

5.2 The pure radiation fluid, k = 0, A = 

In this limit we set w = = 1/3 and 11 = F = 0so that = We conclude 
from p oc that a oc r]. For radiation, the u-equation (75) becomes 

u+{\k^-^)u = 0, (102) 
o r] 



with general solution 



, , . /sin(a;) , A „ / cos(a;) . , ^ \ 
u{x) = A — ^ - cos{x) ] + B { — - sm{x) , (103) 



where we have set x = kr]/^/3 = Cskrj. For the Bardeen potential we obtain 
with (74), up to constant factors, 

^(x) = ^ . (104) 
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On super- horizon scales, a; <C 1, we have 



r/ X ^ 



(105) 



Wc assume that the perturbations have been initiahzed at some early time 
Xin 1 and that at this time the two modes have been comparable. If this is 
the case then B <^ A and we may neglect the -B-mode at later times. 

To determine the density and velocity perturbations and for illustration, 
we also solve the radiation equations using the conservation and Poisson equa- 
tions like for dust. In the radiation case the perturbation equations become 
(with the same notation as above, x = Cgkri) 



4 " 



(106) 

(107) 
(108) 



The general solution of this system is 



cos(x) sin(a;) 

X 



V5 



sin(x) H — cos(a;) 

X 



4-F2a;2 



(109) 
(110) 
(111) 



Again, regularity at a; = requires D\ = 0. Comparing with Eqs. (103,104) 
gives D2 = 2 A. In the super- horizon regime, a; <C 1, we obtain 



3' 



-2A 



3^3 



V = 



A 

— 1=^ 
2V3 



(112) 



On sub-horizon scales, a; 1, we find oscillating solutions with constant 
amplitude and with frequency of k/\/3: 



V 



VSA 



sin(a:) 



Dg = 2Acos(a;) , = -^cos(a;)/a;^ 



(113) 
(114) 



Note that also for radiation perturbations 

A 



D 



3\/3 



x'' 
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is small on super horizon scales, x <C 1. The perturbation amplitude is given by 
the largest gauge invariant perturbation variable. We conclude therefore that 
perturbations outside the Hubble horizon are frozen to first order. Once they 
enter the horizon they start to collapse, but pressure resists the gravitational 
force and the radiation fluid fluctuations oscillate at constant amplitude. The 
perturbations of the gravitational potential oscillate and decay like 1 jc? inside 
the horizon. 

5.3 Adiabatic initial conditions 

Adiabaticity requires that the perturbations of all contributions to the energy 
density are initially in thermal equilibrium. This fixes the ratio of the density 
perturbations of different components. There is no entropy flux and thus F = 
0. Here we consider as a simple example non relativistic matter and radiation 
perturbations. Since the matter and radiation perturbations behave in the 
same way on super-horizon scales, 



D'^J^=A + Bx^, + DC y^"*) oc a;, (115) 



we may require a constant ratio between matter and radiation perturbations. 
As we have seen in the previous section, inside the horizon (x > 1) radia- 
tion perturbations start to oscillate while matter perturbations keep following 
a power law. On sub-horizon scales a constant ratio can thus no longer be 
maintained. There are two interesting possibilities: adiabatic and isocurvature 
perturbations. Here we concentrate on adiabatic perturbations which seem to 
dominate the observed CMB anisotropics. 

Prom r = one easily derives that two components with Pi/pi = 
Wi ^constant, i = 1,2, are adiabatically coupled if (l+wi)!)^^-' = {l+W2)Dg^\ 
Energy conservation then implies that also their velocity fields agree, y^^^ = 
y . This result is also a consequence of the Boltzmann equation in the strong 
coupling regime. We therefore require 



so that the energy flux in the two fluids is coupled initially. 

We restrict ourselves to a matter dominated backgroims, the situation 
relevant in the observed universe after equality. We first have to determine the 
radiation perturbations during a matter dominated era. Since ^ is dominated 
by the matter contribution (it is proportional to the background density of 
a given component), we have \l/ ~ const = !Z'o- We neglect the contribution 
from the sub-dominant radiation to ^. Energy-momentum conservation for 
radiation then gives, with x = kr], 




(116) 



s 3 
y('-)' = 2!f + 



(117) 



(118) 
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Now ^ is just a constant given by the matter perturbations, and it acts 
like a constant source term. The general solution of this system is then 

£)M = Acos{csx) - —Bs\n{c,x) + M [cos(c^x) - 1] (119) 
v3 

y^'^) = Bcos{csx) + ^^sin(csa;) + 2\/3!2>-sin(csa;), (120) 

where Cs = 1/a/3 is the sound speed of radiation. Our adiabatic initial condi- 
tions require 

yir) Y{m) 

lim =Vo= lim < oo. (121) 

x^O X x— >0 X 

Therefore B = and = - 2!Z'. Using in addition ^ = 3Vo (see (101)) 
we obtain 

DW = llfcc.s('^') -8* (122) 

= -^(f> + -x'^) (124) 
" 6 

= iii^x. (125) 

On super-horizon scales, a; <C 1 we have 

20 1 
£)('-) ^ — 1^ and ~ -x^ , (126) 

3 3 

note that D^g^ = (4/3)L»^™^ and V^'') = V^"") for adiabatic initial conditions. 

Another possibility for the initial condition would be iso-curvature initial 
conditions, where you have non- vanishing D'^'^\ /jC"*) and V^^\ y^"^) which 
compensate each other in such a way that if' = on super- horizon scales. 
The simplest inflationary models do not lead to such perturbations and the 
observations imply that they arc not dominating the observed anisotropics 
in the CMB even though they may contribute which could seriously hamper 
the determination of cosmological parameters with CMB anisotropics (see e.g. 
[8, 9] 



6 Scaleir field cosmology 

We now consider the special case of a Priedmann universe filled with self 
interacting scalar field matter. The action is given by 

S = ^ j d'x^\R + j d^x^\ Qa^^^aV - W{^)^ (127) 
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where if denotes the scalar field and W is the potential. The energy momentum 
tensor is obtained by varying the action wrt the metric g^^, 



T,, 



9nu 



The energy density p and the energy flux u are defined by 
For the Friedmann background this gives 



The pressure is given by 



K) = -(i,o) 



1 



We now define the scalar field perturbation, 

ip = (f + d(p . 



(128) 



(129) 



(130) 



(131) 



(132) 



Clearly, the scalar field only generates scalar perturbations (to first order). 
Inserting Eq. (132) in the definition of the energy velocity perturbation v, 



and the energy density perturbation Sp, 

p = p + Sp , 

we obtain 

1 ^ 1 ^2 
6p = -^(fiSip a + W,^ 5ip 



and 



v = - (bix>-\- (3k~'^B) . 



From the stress tensor, Tjj = (yj,j ^p,j — \^^d\^pd^ip + VF] gij we find 

1 ^ 1 ^2 

PT^L = —^^p5w A — W,u, S(p and 11 = . 



(133) 

(134) 
(135) 
(136) 

(137) 



We now define the gauge invariant scalar field perturbation as the value 
of (5<^ in longitudinal gauge. 



(138) 
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The last expression gives Sip^^'^^ in a generic gauge. It is clear that this com- 
bination is gauge -invariant. This variable is very simply related to the other 
gauge-invariant scalar variables. Short calculations give 

V = M(^(s'V'^ (139) 



+ -k-^V 



a 



(140) 



D, = Dg + 3{l + wp (141) 

f>pDs - pSip^s'A (142) 



2W, 



r = 

PP 

n = . (143) 

The last equation shows that the two Bardeen potentials are equal for scalar 
field perturbations, ^ = ^. With this we can write the perturbed Einstein 
equations fully in terms of the Bardeen potential W and V. Since we will 
need them mainly to discuss inflation where curvature plays a minor role, we 
write them down here only for the case of vanishing spatial curvature. Prom 
Eqs. (65) and (67) one can easily generalize to the case with curvature. 

-3H^ -{H + n^- k^)^ = AnGifk-'^iipV + a^W,^ V) (144) 
ip + n^ = AnGip'^k-'^ V (145) 
^ + +{n + 2n^)1r = AnGipk-^i^pV - (pV - c?W,^ V) (146) 

These lead to the following second order equation for the Bardeen potential 
which we have discussed above: 

# + 3(W - (pl(p)if + (2W - m(pl(p + k^)^ = (147) 

or, using the definition = p/p, 

+ m{l + cl)^ + (27i: + (1 + Zcfjl-L^ + k"^)^ = . (148) 

As already mentioned above, this equation differs from the ]!/ equation for a 
perfect fluid only in the last term proportional to fc^. This comes from the 
fact that the scalar field is not in a thermal state with fixed entropy, but it 
is in a fully coherent state {F ^ 0) and field fluctuations propagate with the 
speed of light. On large scales, kt] <^ 1 this difference is not relevant, but on 
sub-horizon scales it does play a certain role. 



7 Generation of perturbations during inflation 

So far we have simply assumed some initial fluctuation amplitude A, with- 
out investigating where it came from or what the fc-dependencc of A might 
be. In this section we discuss the most common idea about the generation of 
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cosmological perturbations, namely their production from the quantum vac- 
uum fluctuations during an inflationary phase. The treatment here is focused 
mainly on getting the correct result with as little eS'ort as possible; we ignore 
several subtleties related, e.g. to the transition from quatum fluctuations of 
the field to classical fluctuations in the energy momentum tensor. The idea is 
of course that the sourse of metric fluctuations are the expectation values of 
the energy momentum tensor operator of the scalar field. 

The basic idea is simple: A time dependent gravitational field very gcner- 
ically leads to particle production, analogously to the electron positron pro- 
duction in a classical, time dependent electromagnetic field. 

7.1 Scalar perturbations 

The main result is the following: During inflation, the produced particles 
induce a perturbed gravitational field with a (nearly) scale invariant spectrum. 



The quantity k^\^{k,r])\^ is the squared amplitude of the metric pertur- 
bation at comoving scale A = 7r/fc. To insure that this quantity is small on a 
broad range of scales, so that neither black holes are formed on small scales nor 
there are large deviation from homogeneity and isotropy on large scales, we 
must require n ~ 1. These arguments have been put forward for the first time 
by Harrison and Zel'dovich [10] (still before the advent of infiation), leading 
to the name 'Harrison- Zel'dovich spectrum' for a scale invariant perturbation 
spectrum. 

To derive the above result we consider a scalar field background dominated 
by a potential, hence a oc |?7|' with g — 1. Looking at the action of this 
system. 



it can be shown (see [11]) that the second order perturbation of this action 
around the Priedmann solution is given by 



k^\^{k,r])\'^ = A;""^ x const, with n ~ 1 . 



(149) 





(150) 



up to some total difi^erential. Here v is the perturbation variable 



a 



C 



(151) 



V = — 



introduced in Eq. (79). Via the Einstein equations, this variable can also 
be interpreted as representing the fluctuations in the scalar field. Therefore, 
we quantize v and assume that initially, on small scales, k\r]\ <C 1, v is in 
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the (Minkowski) quantum vacuum state of a massless scalar field with mode 
function 

Vin = exp{ikr]) . (152) 

The pre-factor vq is a fc-indcpendent constant which depends on convention, 
but is of order unity. From (78) we can derive 

{v/zy = ^ , 

z 

where z cx a is defined in Eq. (81) and u cx arj'l' is given in Eq. (74). On small 
scales, fc|?7| <C 1, this results in the initial condition for u 

Win = exp(ifcr/) . (153) 

The evolution equation for u, (102), reduces in the case of power law 
expansion, a oc |ry|^ to 

9(i±l))„ = 0. (154) 

The solutions to this equation are of the form {k\ri\y/-^ Hj:^\kr]), where /x = 
q+1/2 and Hj^^ is the Hankel function of the ith. kind (i = 1 or 2) of order fj,. 
The initial condition (153) requires that only Hj^^ appears, so that we obtain 



A;3/2 

where again a is a constant of order unity. We define the value of the Hubble 
constant during inflation, which is nearly constant by Hi . With i? = ?i/a ~ 
l/(l^|o) we then obtain a ~ 1/Hi\r]\. Eq. (74) with the Planck mass defiend 



by SttG = Mp^ then gives 



2Mp Mp y I n/ fi \ u \ J 

On small scales this is a simple oscillating function while on large scales k\ri\ <C 
1 it can be approximated by a power law, 

Here we have used /i = l/2 + g<0 and g ~ — 1. This yields 



Mp 

hence n = 1. Detailed studies have shown that even though the amplitude of 'P 
can still be severely affected by the transition from inflation to the subsequent 
radiation era, the obtained spectrum is very stable. Simple deviations from 
de Sitter inflation (like e.g. power law inflation), q > —1 lead to slightly blue 

spectra, n ~ 1. 
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7.2 Vector perturbations 

In the simplest models of inflation where the only degrees of freedom are the 
scalar field and the metric, no vector perturbations are generated. But even 
if they are, subsequent evolution after inflation will lead to their decay. In a 
perfect fluid background, TJy = 0, vector perturbations evolve according to 
Eq.(71) which implies 

r^oca^'^'-^ (158) 

For a radiation fluid, p/ p = c?^ < 1/3, this leads to a non-growing vorticity. 
The dynamical Einstein equation (68) gives 

(159) 

and the constraint (66) reads (at early times, so that we can neglect curvature) 

n^{kr]fa^^l (160) 

Therefore, even if they are created in the very early universe on super-horizon 
scales during an inflationary period, vector perturbations of the metric decay 
and become soon entirely negligible. Even if Hi remains constant in a radiation 
dominated universe, it has to be so small on relevant scales at formation 
{kr]in "C 1) that we may safely neglect it. 

Vector perturbations are irrelevant if perturbations have been created at 
some early time, e.g. during inflation. This result changes completely when 
considering 'active perturbations' like for example topological defects where 
vector perturbations contribute significantly to the CMB anisotropics on large 
scales, see Ref. [12]. Furthermore, it is interesting to note that vector pertur- 
bations do not satisfy a wave equation and therefore will in no case show 
oscillations. Vorticity simply decays with time. 



7.3 Tensor perturbations 

The situation is different for tensor perturbations. Again we consider the 
perfect fluid case, n^J^ = 0. Eq. (69) implies, if k is negligible, 

Hij + ^Hij + k^Hij = . (161) 

If the background has a power law evolution, a rf^ this equation can be 
solved in terms of Bessel or Hankel functions. The less decaying mode solution 

to Eq. (161) is Hij — eijX^/"^"^ Ji/2-q{x), where Ji, denotes the Bessel function 
of order v, x = krj and e^- is a transverse traceless polarization tensor. This 
leads to 

Hij = const for a; <§; 1 (162) 
Hij = - for a; ~ 1 . (163) 
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One may also quantize the tensor fluctuations which represent gravitons. Do- 
ing this, one obtains (up to log corrections a scale invariant spectrum of tensor 
fluctuations from inflation: For tensor perturbations the canonical variable is 
simple given by hij = MpaHij. The evolution equation for hij = heij is of 
the form 

'h+{P +m'^{r]))h = , (164) 

where w?{r]) = —a/ a. During inflation = —q{q— 1) is negative, leading to 
particle creation. Like for scalar perturbations, the vacuum initial conditions 
are given on scales which are inside the horizon, k'^ ;:g> |m^|, 

hin — —j=exp{kr]) for fc|77| » 1. 
Vk 

Solving Eq. (164) with this initial condition, gives 

/^=-^(fc|r;|)V2i^(^\/^(fcr,), 

where Hi^^ is the Hankel function of degree u of the second kind. On super 
horizon scales, H^^^^^{kr]) cx (fc|?7|)'^~-'-/^ this results in ~ |?7|(fc|?7|)^'^~-'-. 
Using the relation between hij = hcij and Hij one obtains the spectrum 
of tensor perturbation generated during inflation. For exponential inflation, 
g ~ — 1 one finds again a scale invariant spectrum for Hij on super-horizon 
scales. 

k^\HijW^c^{Hi^/Mpf ock""^ with nr . (165) 



8 Lightlike geodesies and CMB anisotropies 

After decoupling, rj > r^deo photons follow to a good approximation light-like 
geodesies. The temperature shift of a Planck distribution of photons is equal 
to the energy shift of any given photon, which is independent of the photon 
energy (gravity is 'achromatic'). 

The unperturbed photon trajectory follows 

rvo 

{x''{'n)) = {'n, n(7?')rfr?'+xo), 

Jri 

where xq is the photon position at time rjo and n is the (parallel transported) 
photon direction. With respect to a geodesic basis (e)^^^, the components of 
n are constant. If k = we may choose = d/dx^; if k 7^ these vector 
fields are no longer parallel transported and therefore do not form a geodesic 
basis (Ve^ej = 0). 

Our metric is of the form 
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ds^ = a^ds^ ,with (166) 
ds'^ = {'yij,^ + h^^)dx^'dx'', 700 = -1, 7^0 = 0, 7^- = 7^, (167) 

as before. 

We make use of the fact that Ught-Uke geodesies are conformally invari- 
ant. More precisely, ds^ and ds"^ have the same hght-hke geodesies, only the 
corresponding affine parameters are different. Let us denote the two affine 
parameters by A and A respectively, and the tangent vectors to the geodesic 
by 

n = n = -n = ^ , = = , n° = 1 , = 1. (168) 

d\ a dX 

We set n° = 1 + 5n°. The geodesic equation for the perturbed metric 

ds^ = [ly^u + h^u)dx''dx'' (169) 

yields, to first order, 

J-^n" = -Jr^^n^n'^. (170) 

For the energy shift, we have to determine S'nP. Since 5°'* = —1 • Sq^ + 
first order, we obtain SF^^ = — V2(/ic«o|/3 + /^/JOla ~ /^a/s), so that 



^feO = /iaO|/3«''ri" - i/ia/3n"n^. (171) 



Integrating this equation we use hao\f3n^ = -jx(^"On"), so that the change of 
nP between some initial time r]i and some final time r]f is given by 

5n°|{ = [/loo + hojTi^] i-\f^ h^^nVdX . (172) 

On the other hand, the ratio of the energy of a photon measured by some 
observer at to the energy emitted at ti is 

Ef _ {n- u)f _Tf {n- u)f 
Ei {n ■ u)i Ti (n ■ u)i 

where u / and Ui are the four- velocities of the observer and emitter respectively, 
and the factor Tf/Ti is the usual (unperturbed) redshift, which relates n and 
n. The velocity field of observer and emitter is given by 

u = (1 - A)dr^ + v'di . (174) 

An observer measuring a temperature To receives photons that were emit- 
ted at the time rjdec of decoupling of matter and radiation, at the fixed tem- 
perature Tdec- In first-order perturbation theory, we find the following relation 
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between the unperturbed temperatures Tf, Tj, the measurable temperatures 
To = Tj + 5Tf, T^ec = Ti + STi, and the photon density perturbation: 



Tf _ To STf . 5T,\ T, 



(175) 



where 5^''^ is the intrinsic density perturbation in the radiation and we have 

used in the last equality. Inserting the above equation and Eq. (172) 

into Eq. (173), and using Eq. (30) for the definition of /i^i^, one finds, after 
integration by parts [6] the following result for scalar perturbations: 



Ef _ Tq 



+ V^n^ + + # \ + S)d)\ . (176) 



Here o'f^ denotes the density perturbation in the radiation fluid, and V'^^^ 
is the peculiar velocity of the baryonic matter component (the emitter and 
observer of radiation). 

Evaluating Eq. (176) at final time ?7o (today) and initial time r^dec, we 
obtain the temperature difference of photons coming from different directions 

f .^-^.|(n,-|(nO. (177) 

Direction independent contributions to do not contribute to this difference. 
We also do not want to include the term V,- (770)71-' which simply describes the 
dipole due to our motion with respect to the emission surface and which is 
much larger than the contributions from the higher multipoles. Therefore we 
can set 



ZiT(n) 



ect ^dec 

(178) 

where x{r]) is the unperturbed photon position at time r] for an observer at 

xo, and x^ec = x{ridec) (If k = we simply have x(77) = xq - (770 - ??)"•)• The 
first term in Eq. (178) is the one we have discussed in the previous chapter. 
It describes the intrinsic inhomogeneities on the surface of last scattering, 
due to acoustic oscillations prior to decoupling. Depending on the initial con- 
ditions, it can contribute significantly also on super-horizon scales. This is 
especially important in the case of adiabatic initial conditions. As we have 
seen in Eq. (126), in a dust -I- radiation imivcrsc with J? = 1, adiabatic initial 
conditions imply D^''\k, rj) = -f^{k, 77) and V'^'>'> = < D^"'^ for krj < 1. 
With ^ = ^ the the square bracket of Eq. (178) therefore gives for adiabatic 
perturbations 

V / adiabatic 
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on super- horizon scales. The contribution to ^ from the last scattering sur- 
face on very large scales is called the 'ordinary Sachs- Wolfe effect' (OSW). 
It has been derived for the first time by Sachs and Wolfe [13] in 1967. For 
isocurvature perturbations, the initial condition olf^ {k, 77) ^ for — > is 
satisfied and the contribution of Dg to the ordinary Sachs- Wolfe effect can be 
neglected. 



isocurvature 

The second term in (178) describes the relative motion of emitter and observer. 
This is the Doppler contribution to the CMB anisotropics. It appears on the 
same angular scales as the acoustic term; we call the sum of the acoustic and 
Doppler contributions "acoustic peaks" . 

The last two terms are due to the inhomogencities in the spacetime geom- 
etry; the first contribution determines the change in the photon energy due 
to the difference of the gravitational potential at the position of emitter and 

(r) 

observer. Together with the part contained in Dg they represent the "or- 
dinary" Sachs- Wolfe effect. The integral accounts for red-shift or blue-shift 
caused by the time dependence of the gravitational field along the path of the 
photon, and represents the so-called integrated Sachs- Wolfe (ISW) effect. In 
a n — I, pure dust universe, the Bardeen potentials are constant and there 
is no integrated Sachs- Wolfe effect; the blue-shift which the photons acquire 
by falling into a gravitational potential is exactly canceled by the redshift 
induced by climbing out of it. This is no longer true in a universe with sub- 
stantial radiation contribution, curvature or a cosmological constant. 

The sum of the ordinary Sachs-Wolfe term and the integral is the full 
Sachs- Wolfe contribution (SW). 

For vector perturbations 5^'^^ and A vanish and Eq. (173) leads to 

(i?//i?,)(^) = ia,/af)[l - y/'-^n^lf + J^^ a.n^dX] . (179) 

We obtain a Doppler term and a gravitational contribution. For tensor per- 
turbations, i.e. gravitational waves, only the gravitational part remains: 

rf 

(Ef/Ei)^'^^ ={ai/af)[l- Hijn^n^dX] . (180) 

J i 

Equations (176), (179) and (180) are the manifestly gauge invariant results 
for the energy shift of photons due to scalar, vector and tensor perturba- 
tions. Disregarding again the dipole contribution due to our proper motion, 
Eqs. (179,180) imply the vector and tensor temperature fluctuations 



AT{n) 



= Vl"'\vdec,Kdec)n' + &j{v,Av))n'dX (181) 
(^^^Y^ =-l^^ Hij{vMv)Wn^dX . (182) 
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Note that for models where initial fluctuations have been laid down in the 
very early universe, vector perturbations arc irrelevant as we have already 
pointed out. In this sense Eq. (181) is here mainly for completeness. However, 
in models where perturbations are sourced by some inherently inhomogeneous 
component {e.g. topological defects, see Ref. [12]) vector perturbation can be 
important. 

9 Power spectra 

One of the basic tools to compare models of large scale structure having 
stochastic initial fluctuations with observations are power spectra. They are 
the "harmonic transforms" of the two point correlation functions. If the per- 
turbations of the model under consideration are Gaussian (a relatively generic 
prediction from inflationary models), then the power spectra contain the full 
statistical information of the model. 

Let us first consider the power spectrum of dark matter, 



Here ( ) indicates a statistical average, ensemble average, over "initial condi- 
tions" in a given model. Pd (k) is usually compared with the observed power 
spectrum of the galaxy distribution. This is clearly problematic since it is by 
no means evident what the ratio of these two spectra should be. This problem 
is known under the name of 'biasing' and it is very often simply assumed that 
the dark matter and galaxy power spectra differ only by a constant factor. 
The hope is also that on suflicently large scales, since the evolution of both, 
galaxies and dark matter is governed by gravity, their power spectra should 
not differ too much. This hope seems to be reasonably justifled [14]. 
The power spectrum of velocity perturbations satifies the relation 



For ~ we have used that |fcy|(77o) = Dg"^\r]o) ~ Hofl'^-^Dg on sub-horizon 
scales (sec e.g. [15]). 

The spectrum we can be both, measured and calculated to the best accu- 
racy is the CMB anisotropy power spectrum. It is defined as follows: AT jT is 
a function of position xq, time t/o and photon direction n. We develop the n- 
dependence in terms of spherical harmonics. We will suppress the argument 
r\Q and often also xq in the following calculations. All results are for today 
(?7o) and here (xo). By statistical homogeneity statistical avc^ra.gc^s over an 
ensemble of realisations (expectation values) are supposed to be independent 
of position. Furthermore, we assume that the process generating the initial 
perturbations is statistically isotropic. Then, the off-diagonal correlators of 
the expansion coefficients vanish and we have 




(183) 




(184) 
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— (xo,n,r7o) = ^ a^m(xo)yem(n), (a^„ • a|,„,) = (5«'(5mm'Q (185) 

The Q's are the CMB power spectrum. 

The two point correlation function is related to the Q's by 

(^(n)^(n')) = X! {o-irn- a},^,)Yi,^{n)Yl,^,{n') = 



^ 1 
53 Q ^ r,„(n)r;„(n') = — ^(21 + i)qp,(m), (186) 

2|±iP.(n.n') 

where we have used the addition theorem of spherical harmonics for the last 

equality; the P^'s are the Legendre polynomials. 

Clearly the ajm's from scalar, vector and tensor perturbations are uncor- 
related, 

a^s) (V) \ _ (AS) (T) \ _ /AV) (T) \ _ Q .^g„. 

Since vector perturbations decay, their contributions, the C|^\ are neg- 
ligible in models where initial perturbations have been laid down very early, 
e.g. , after an inflationary period. Tensor perturbations are constant on super- 
horizon scales and perform damped oscillations once they enter the horizon. 

Let us first discuss in somewhat more detail scalar perturbations. We spe- 
cialize to the case k = for simplicity. We suppose the initial perturbations 
to be given by a spectrum, 

We multiply by the constant 'r]'^~^ , the actual comoving size of the horizon, 
in order to keep A dimensionless for all values of n. A then represents the 
amplitude of metric perturbations at horizon scale today, k = l/rjo- 
On super-horizon scales we have, for adiabatic perturbations: 

= + 0{{kr)f), V^''^ = F^'-) = 0(kr]) (189) 

4 3 

The dominant contribution on super-horizon scales (neglecting the inte- 
grated Sachs-Wolfe effect / ^ — ) is then 

AT 1 

— (xo,n,7?o) = ^'^(^dec,??dec)- (190) 

The Fourier transform of (190) gives 
AT 1 

— (k, n,r,o) = -^{k, %ec) • e*"^'"'-''^-) . (191) 
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Using the decomposition 

oo 

gikn(,,o-r,aec) = Y^^2i + l)i'je{k{m - %ec))PKk " n) , 

where jt are the spherical Bessel functions, we obtain 
AT AT \ 

— (xo,n,r?o) — (xo,n',r?o)^ (192) 



If I AT AT 
= — / d^a;o ^— (xo,n,r7o)^(xo,n',r;o) 

1 f . / AT f AT\* 

" (2^ J — (k,n,r?o) J (k,n',r?o) 

^ ' •' i,e'=o 

■je{k{vo - Vdcc))je'{k{r]o - ??doc))-P£(k • n) • Pe,{k ■ n') . (193) 

In the second equal sign wo have used the unitarity of the Fourier transfor- 
mation. Inserting Peikn) = Y.m ^/m(k)^«m(n) and 
P^'(kn') = 2^^!^™' ^£*m'(k)^<?'m'(n'); integration over the directions dOj, 

gives Su'Smm' J2m. ^/m ("i^^™ ) • 

Using as well Y.m ^^m(")^^™("') = where /i = n • n', we find 

/AT AT , \ 

( ^(xo,n,r/o)^(xo,n,?7o) > 

^ ^P.(m)^ Jf{l\^\') k'mvo - mec)). (194) 

Comparing this equation with Eq. (186) we obtain for adiabatic perturba- 
tions on scales 2 < £ <^ xivo — ??dec)/%ec ~ 100 



n(SW) ^ ^(OSW) 



2 r dk ' ' 



k 



k'n (k ivo - Vdec)) . (195) 



If is a pure power law as in Eq. (188) and we set k(j]Q — r/jec) ~ krjo, the 
integral (195) can be performed analytically. For the ansatz (188) one finds 



9 23-"r2(2- f)r(f + 1 - f) 



Ci^w) ^ ^" 2 '2> for - 3 < n < 3 . (196) 



Of special interest is the scale invariant or Harrison-Zel'dovich spectrum, 
n = 1 (see Section 7). It leads to 
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(SW) _ 



= const. ~ 



= IT/ 



(197) 



This is precisely (within the accuracy of the experiment) the behavior observed 
by the DMR experiment aboard the satellite COBE [16]. 

Inflationary models predict very generically a HZ spectrum (up to small 
corrections). The DMR discovery has therefore been regarded as a great suc- 
cess, if not a proof, of inflation. There are other models like topological defects 
[17, 18, 19] or certain string cosmology models [20] which also predict scale- 
invariant, i.e. , Harrison Zel'dovich spectra of fluctuations. These models do 
however not belong to the class investigated here, sinc;c; in these models per- 
turbations are induced by seeds which evolve non-linearly in time. 

For isocurvature perturbations, the main contribution on large scales 
comes from the integrated Sachs-Wolfe effect and (195) is replaced by 



(ISW) 



(198) 



Inside the horizon is roughly constant (matter dominated). Using the ansatz 
(188) for ^ inside the horizon and setting the integral in (198) ^ 2'P{k, rj = 
l/k)j^{kr]o), we obtain again (196), but with replaced by 4A^. For a 

fixed amplitude A of perturbations, the Sachs- Wolfe temperature anisotropics 
coming from isocurvature perturbations are therefore about a factor of 6 times 
larger than those coming from adiabatic perturbations. 

On smaller scales, £ ~ 100 the contribution to AT/T is usually dominated 
by acoustic oscillations, the first two terms in Eq. (178). Instead of (198) we 
then obtain 



{AC) 



2 r°° dk / 1 



. (199) 



To remove the SW contribution from Dg ' we have simply replaced it 
by D'^^^ which is much smaller than ^ on super-horizon scales and therefore 



does not contribute to the SW terms. On sub-horizon scales D'^^'' ~ D 



(r) 
9 

and y^'"^ are oscillating like sine or cosine waves depending on the initial 
conditions. Correspondingly the cj"^'^'' will show peaks and minima. On very 
small scales they are damped by the photon diffusion which takes place during 
the recombination process (see contribution by A. Challinor). 

For gravitational waves (tensor fluctuations), a formula analogous to (196) 
can be derived, 



dp = ^ J dkk"^ / J"" dr)H{r], k) 



JdHvo - v)) 



2)! 



(£-2)!- 



(200) 
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Fig. 1. A COBE normalized sample adiabatic (solid line) and isocurvature (dashed 

line) CMB anisotropy spectrum, £{£+l)Ci, are shown on the top panel. The quantity 
shown in the bottom panel is the ratio of temperature fluctuations for fixed value 
of A (from Kanazawa et al. [21]). 



To a very crude approximation we may assume H = on super-horizon 
scales and / dr]Hji{k{riQ — rj)) ~ H{r] = l/k)ji{krio). For a pure power law, 

{\H{k,v = l/kf) = ^|A;"-r?o (201) 

one obtains 

,T)^2i£ + 2l fdx jfix) 



(f + 2)!^2 r{<o-nT)r{l~2 



2 ' 



(202) 



{e-2)\ ^2*^-^-rn{l-nT)r{£ + A-^)' 
For a scale invariant spectrum {ht ~ 0) this results in 

The singularity at ^ = 2 in this crude approximation is not real, but there is 
some enhancement of £{£ + 1)C^^^ at ^ 2 see Fig. 2). 

Since tensor perturbations decay on sub-horizon scales, i ^ 60, they are 
not very sensitive to cosmological parameters. 
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Again, inflationary models (and topological defects) predict a scale invari- 
ant spectrum of tensor fluctuations {n^ ^ 0). 

On very small angular scales, £ ~ 800, fluctuations are damped by col- 
lisional damping (Silk damping). This effect has to be discussed with the 
Boltzmann equation for photons which is presented in detail in the course by 
A. ChaUinor. 

Scalar Tensor 




10 100 IQOO 10 100 IQOO 



Fig. 2. Adiabatic scalar and tensor CMB anisotropy spectra are shown (top panels). 
The bottom panels show the corresponding polarization spectra (see course by A. 
Challinor). (from [22]). 



10 Some remarks on perturbation theory in braneworlds 

Since there has been so much interest in them recently, let me flnally make 
some remarks on perturbation theory of five dimensional braneworlds. I shall 
just present some relatively simple aspects without derivation. A thorough 
discussion of braneworlds is given in the course by R. Maartens. Different 
aspects of the perturbation theory of braneworlds can be found in the growing 
literature on the subject [23, 24]. 

The bulk background metric of a flve dimensional braneworld has three 
dimensional spatial slices which homogeneous and isotropic, hence spaces of 
constant curvature. Its line element is therefore of the form 
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ds^ = -n^(t, y)df + a^it, y)-f,jdx'dx^ + b^it, y)dy^ = gAsdx^dx^ (204) 

Perturbations of such a spacetime can be decomposed into scalar, vector and 
tensor modes with respect to the three dimensional spatial slices of constant 
curvature. One can always choose the so-called generalized longitudinal 
gauge such that the perturbations of the metric are given as follows: 

ds^ = -71^(1 + 2^)dt^ + [(1 - 2<Z')7„- + 2Hij] dx'dx^ + b^{l + 2C)dy^ 

-2nbBdtdy - 2naSidx'dt + ab2£,dx'dy . (205) 

Here Si and are divergence free vector fields, vector perturbations and 
Hij is a divergence free traceless symmetric tensor, the tensor perturbation. 
'F. <P, C and B arc four scalar perturbations. 

It can be shown that this choice determines the gauge completely. One 
can actually define gauge invariant perturbation variables which reduce to 
the ones above in the generalized longitudinal gauge [24]. Writing down the 
perturbed Einstein equations for these variables in the most general case is 
quite involved. These equations can be found in Ref. [24], but I don't want 
to repeat them here. I just discuss their general structure in the ease of an 
empty bulk. It is clear that \l/ and ^ correspond to the Bardeen potentials of 
four dimensional cosmology. Si is the four dimensional vector perturbation 
and Hij represents four dimensional gravitational waves. C and B as well as 
£i are new degrees of freedom which are not present in the four dimensional 
theory. 

If we assume vanishing perturbations of the bulk energy momentum tensor 
e.g. if the bulk is anti de Sitter like in the Randall Sundrum model [25] called 
RSII in what follows (see course by R. Maartens), the perturbation equations 
reduce to 

□5(!f- + ^) = (206) 
a^Si = (207) 
OsHij = (208) 

and all the other perturbation variables are determined by constraint equa- 
tions. Here Ds is the five dimensional d'Alembertian with respect to the back- 
ground metric (204) . This structure of the equations is to be expected: Gravi- 
tational waves in d spacetime dimensions are a spin 2 field with respect to the 
group of rotations SO{d — 2) since they are massless (see e.g. [26]). For d = 5, 
d—2 = 3 they therefore have 5 degrees of freedom. These correspond exactly 
to the one scalar (206), two vector (207) and two tensor (208) degrees of free- 
dom with respect to the 3-dimensional slices of conrant curvature. These free 
massless degrees of freedom obey the wave equations above. 
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The perturbed Israel junction conditions (see course by R. Maartens) then 
determine boundary conditions for the behavior of the perturbations at the 
brane position(s). 

As an example, I write down the vector and tensor perturbation equations 
and their bulk solutions for the RSII model which has only one branc. There 
the bulk is a five dimensional anti-de Sitter spacetime, and we can choose 
coordinates so that the background metric has the form 



ds^ = [-dt^ + Sijdx'dx^ + dy^ 



(209) 



and 

□5 = -d1 + A + dl-^dy , (210) 

where A denotes the three dimensional spatial Laplacian. For an arbitrary 
mode which satisfies this wave equation we make the ansatz 

(t>{t, X, y) = exp(i(k • x - ujt))(j){uj, k, y) . 

We then obtain a Bessel differential equation for <p{co, k, y) with general solu- 
tion 

(j) = j4(a;,k)(ym)^J2(my) + B(a;,k)(ym)^l2(my) , = u"^ — k'^ , 

= (l>A + cl>B. (211) 

These modes are normalizable in the sense that 



/•CO 

/ \<l>\'^V-9dy < 00 . 
Jvb 



Here i/6 > is the branc position. The Israel junction condition for the vector 
and tensor modes for RSII become 

-dyHij {y = yb) = nlHij (yt) (212) 
-Eiiyb) = 4ni{yb) , (213) 

where ily respectively TTj are the tensor rsp. vector contribution to the 
anisotropic stresses on the brane. If the latter vanish, the junction condition 
simple requires 

B{co,k) = A{cj,k)J2imyb)/Y2imyb) . 

This result has been derived for the tensor mode by Randall and Sundrum [25] . 
For scalar perturbations the situation is somewhat more complicated since 
there is an additional degree of freedom which is the perturbed position of 
the brane, yt ^ yb + e- It is rather subtle to take this brane bending correctly 
into account. A very interesting work showing that this effect is actually most 
relevant to obtain the correct Newtonian limit in the RSII model can be found 
in Ref. [27] 
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Let us also briefly discuss the zero-mode. Prom the brane point of view, 
the modes diseussed here represent waves (particles) which couple only to the 
energy momentum tensor of the brane and which obey a dispersion relation 

— fc^ = m^, hence the parameter m of the solutions (213) is their mass. In 
the limit m — > the solutions turn into power laws in y, 



Of these modes, for tensor and vector perturbations, the B mode is normaliz- 
able. For scalar perturbations the situation is more complicated. For a mode 
to be 'normalizable' we want all the perturbations to be normalizable. 



But from the constraints one obtains for m = C oc which is not nor- 
malizable. This mode diverges logharitmically and only converges due to the 
oscillations of the Bessel functions if m 7^ 0. 

If all the five graviton zero-modes are normalizable, like in all compact 
braneworlds, e.g. in models with two branes, the vector and scalar mode lead 
to three additional degrees of freedom (to the usual two four-dimensional 
graviton modes) which couple via the junction conditions to the brane energy 
momentum tensor and spoil the phenomenology of the model. The resulting 
four dimensional gravity is not Einstein gravity but Brans-Dicke or even more 
complicated. As an example, in Ref. [28] the contribution from the scalar zero- 
mode (206) to the gravity wave emission from a binary pulsar is calculated 
and shown to be in contradiction with observations. This problem is well- 
known from Kaluza Klein theories, it is the so called moduli problem. The 
way out of it is usually to render the modes massive. There are different 
suggestions how this can be achieved for the scalar mode. One possibility is 
the so called Goldberger-Wise mechanism [29] which shows that under certain 
circumstances a bulk scalar field can do the job. But certainly, a physically 
acceptable braneworld is defined only together with its mechanism how to get 
rid of such unwanted modes (see also Ref. [30] . 

The advantage of the RSII model is that the scalar gravity wave zero- 
mode is not normalizable in this model and therefore it does not contribute. 
This very promising property of the RSII model has let to its popularity. It 
seems to induce the correct four dimensional Einstein gravity on the brane, 
in the cosmological context this property is still maintained at sufficiently low 
energies. 

11 Conclusions 



My) = Ay^ + B . 



(214) 




Vb 



for = ^ , W , C and B . 



In this course I have; given an introduction to cosmological perturbation the- 
ory. Perturbation theory is an important tool especially to calculate CMB 
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anisotropics and polarisation since these are very small and can be deter- 
mined reliably within linear cosmological perturbation theory. To determine 
the evolution of the cosmic matter density, linear perturbation theory has to 
be complemented with the theory of weakly non-linear Newtonian gravity and 
with N-body simulations. To finally understand the formation of galaxies non- 
gravitational highly non-linear physics, like heating and cooling mechanisms, 
dissipation, nuclear reactions etc. have to be taken into account. This very 
difficult subject is still in its infancy. 

To make progress in our understanding of braneworlds, linear perturbation 
theory can also be most helpful. We can use it to determine e.g. the propagat- 
ing modes of the gravitational field on the brane, light deflection and redshift 
in weak gravitational fields and the Newtonian limit. The condition that linear 
perturbations on the brane at low energy and large distances reduce to those 
resulting from Einstein gravity is non-trivial and has, to my knowledge, not 
yet been fully explored to limit braneworld models. 

Acknowledgement. I thank the organizers for a well structured school in a most 
beautiful environment. 
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